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We present a measurement of the time-dependent CP-violating asymmetries in B° — > K^ir 
decays based on 383 million T(45) — > BB events collected by the BABAR, experiment at the PEP-II 
asymmetric-energy B Factory at SLAC. We measure the direct CP-violating asymmetry C K o T o = 
0.24 ± 0.15 ± 0.03 and the CP- violating asymmetry in the interference between mixing and decay 
SjfO w o = 0.40 ± 0.23 ± 0.03, where the first errors are statistical and the second are systematic. On 
the same sample, we measure the decay branching fraction, obtaining B(B° — » K°ir°) = (10.3 ± 
0.7 ±0.6) x 10 -6 . 

PACS numbers: 13.25.Hw, 11.30.Er 



The BABAR and Belle experiments have measured the 
weak phase /3 [H, Q of the Cabibbo-Kobayashi-Maskawa 
(CKM) quark mixing matrix Q with a better preci- 
sion than the Standard Model (SM) prediction [4( de- 
rived from measurements of other CP-conserving and 
CP-violating processes. The agreement between the the- 
oretical and experimental results has shown that the 
CKM matrix correctly describes these measurements of 
(3 to good precision. 

A major goal of the B Factory experiments is now to 
search for indirect evidence of New Physics (NP). One 
strategy is to compare the measured value of the CP 
violation (CPV) parameters from b — > sec to indepen- 
dent determinations of the same quantities using pro- 
cesses that are sensitive to the contributions of NP effects 
through loop diagrams. 

CPV in B decays to a final state / can be parame- 
terized by Cf, measuring direct CPV, and Sf, measur- 
ing CPV in the interference between decays with and 
without mixing. In the Standard Model for penguin- 
dominated processes b — > sqq (q — u,d,s) [2lj, Sf and 
Cf are expected to be consistent with the values from 
b — ► sg decays. Additional CKM suppressed contribu- 
tions to the amplitude can induce only small deviations 
from this expectation. On the other hand, additional 
loop contributions from NP processes may produce ob- 
servable deviations [f| Q . 

The CKM and color suppression of the tree-level b — ► 
suu transition leads to the expectation that the de- 
cay B° — > K°tt° is dominated by a top quark medi- 
ated b — > sdd penguin diagram, which carries a weak 
phase arg(V(bV t *). If non-leading contributions are small, 
S K o^ is expected to be equal to sin 2/3 and C K a v o ~ 0. 

In addition, it is possible to combine the direct CP 
asymmetries and the branching fractions of the four 
B — > Kit modes to test precise sum rules @, H, Hj ■ The 
experimental uncertainty on these sum rules is dominated 
by the error on the direct CP asymmetry in B° — > K°ir°. 
Therefore a precise measurement of both the direct CP 
asymmetry, and the branching fraction, in this decay 
channel represents an important consistency test of the 
SM. 

The time-dependent CP asymmetries of the decay 
B° — ► K®tt (K® — > 7r + 7r~) have been measured by 
BABAR 10] and subsequently by Belle 11 1, and both ex- 
periments have also measured the branching ratio 

0, 



Il3| . In this work, we present an update of the these re- 
sults based on 383 million Y(4S) — > BB decays collected 
with the BABAR detector at the PEP-II e + e~ collider, 
located at the Stanford Linear Accelerator Center. 

The BABAR detector, which is described elsewhere flij ]. 
provides charged particle tracking through a combination 
of a five-layer double-sided silicon micro-strip detector 
(SVT) and a 40-layer central drift chamber, both op- 
erating in a 1.5 T magnetic field to provide momentum 
measurements. Charged kaon and pion identification is 
achieved through measurements of particle energy loss 
in the tracking system and Cherenkov cone angle in a 
detector of internally reflected Cherenkov light. A seg- 
mented CsI(Tl) electromagnetic calorimeter (EMC) pro- 
vides photon detection and electron identification. Fi- 
nally, the instrumented flux return of the magnet allows 
discrimination between muons and pions. 

We reconstruct if' — > 7r + 7r~ candidates from pairs of 
oppositely charged tracks. The two-track combinations 
must form a vertex with a \ 2 probability greater than 
0.001 and a ir + ir~ invariant mass within 11.2 MeV/c 2 
(3.7a) of the Kg mass [l5|. We form 7r° — > 77 candi- 
dates from pairs of energy depositions in the EMC that 
are isolated from any charged tracks, carry a minimum 
energy of 50 MeV per photon, fall within the mass win- 
dow 110 < to 77 < 160 MeV/c 2 , and have the expected 
lateral shower shapes. Finally, we construct B° — > K®ir° 
candidates by combining K® and ir° candidates in the 
event using kinematic and geometric information of the 
decay which constraints the B° decay vertex to originate 
in the e + e~ interaction region. We extract the flight 
length of the K® from the fit and require that the re- 
constructed proper lifetime be greater than five times its 
uncertainty. We require that the x 2 probability of the fit 
be greater than 0.001. 

For each B° candidate two, independent kinematic 
variables are computed. The first one is ms, the in- 
variant mass of the reconstructed B meson, Bcp- The 
second one is m m ; ss , the invariant mass of the other B, 
^tagj computed from the known beam energy, apply- 
ing a mass constraint to Bcp [Hj]. For signal decays, 
n%B (fimiss) peaks near the B° mass with a resolution 
of - 36MeV/c 2 (~ 5.3 MeV/c 2 ). Both the m miss and 
mj distributions exhibit a low-side tail from leakage of 
energy deposits out of the EMC. We select candidates 
within the window 5.11 < m m i ss < 5.31 GeV/c 2 and 
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5.13 < rns < 5.43 GeV/c 2 , which includes the signal 
peak and a "sideband" region for background characteri- 
zation. For the 0.8 % of events with more than one recon- 
structed candidate, we select the combination with the 
smallest x 2 = X^",*: ( TO * ~ TO ';) 2 / cr m I > where rrn (m^) 
is the measured (nominal) mass and a mi is the estimated 
uncertainty on the measured mass of particle i. 

We exploit topological observables, computed in the 
T(4S) rest frame, to discriminate jet-like e + e~ to 
qq events (q — u, d, s, c) from more spherical BB 
events. We compute the value of L^JLq, where Lj = 
J2i I Pi II cos 6** p . Here, p* is the momentum of parti- 
cle i and 8* is the angle between p* and the sphericity 
axis [III of the Bcp candidate, and the sum does not in- 
clude the decay tree of the -Btag. In order to reduce the 
number of background events, we require L 2 /Lq < 0.55. 
We compute cos 9* B , the cosine of the angle between the 
direction of the B meson and the nominal direction of 
the magnetic field (z axis). This variable is distributed 
as 1 — cos 2 9* B for signal events and nearly flat for back- 
ground events. We select events with | cos 6*^ | < 0.9. 
We also use the distribution of L 2 /L and of cos#^ to 
discriminate the signal from the residual background in 
a maximum likelihood fit. Using a full detector simula- 
tion, we estimate that our selection retains (33.6 ± 1.6)% 
of the signal events, where this error includes statisti- 
cal and systematic contributions. The selected sample of 
B° — > K®tt° candidates is dominated by random K®tt° 
combinations from e + e~ — > qq (q — u, d, s, c) fragmen- 
tation. Using large samples of simulated BB events, we 
find that backgrounds from other B meson decays can 
be generally neglected, but we include some specific B 
decay channels in our study of the systematic errors. 

For each B° — > K®ir° candidate, we examine the re- 
maining tracks and neutral candidates in the event to 
determine if the -B tag meson decayed as a B° or a B° 
(flavor tag). We use a neural network to determine the 
flavor of the Stag meson from kinematic and particle 



the distribution of the difference of the proper decay 
times, At = tcp — ftagi where the tcp refers to the Bcp 
and t tag to the Stag- At the T(AS) resonance, the At 
distribution follows 



identification information [17| . Each event is assigned 



to one of six mutually exclusive tagging categories, de- 
signed to combine flavor tags with similar performance 
and vertex resolution. We measure the performance of 
this algorithm in a data sample (Bfl av ) of fully recon- 
structed B° — > D(*'7r + /p + /di decays. The average 
effective tagging efficiency obtained from this sample is 

Q = E c e s(! - 2wC ) 2 = ( 30 - 5 ± °- 3 ) % > where e s and wC 

are the signal efficiency and mistag probability, respec- 
tively, for events tagged in category c, and the error is 
statistical only. We take into account differences in tag- 
ging efficiency (for signal and background) and mistag 
(only for signal) for B° and B° events, in order to exclude 
any source of fake CPV effects. For the background, the 
fraction of events (e c B ) and the asymmetry in the rate of 
B" versus B° tags in each tagging category are extracted 
from the fit to the £?fl av data. 

Time-dependent CP asymmetries are determined from 



V{At) = 



,-\At\/r 

At 



x{1 ± [ S f sin (Am d At) - C f cos (Am d At) } } , (1) 

where the +(-) sign corresponds to B tag decaying as B° 
(B°), t is the neutral B lifetime, Arrid is the mixing angu- 
lar frequency, C/ is the magnitude of direct CP violation 
in the decay to final state /, and Sf is the magnitude 
of CP violation in the interference between mixing and 
decay. To account for flavor mistags we reduce Sf by the 
factor 1 — 2w c . For the case of penguin dominance, we 
expect Sjfo^o ~ sin 2/3, and Cr-o^o ~ 0. 

The reconstructed proper time difference At r is com- 
puted from the measured Az = zep — ztagj the difference 
of the reconstructed decay vertex positions of the Bcp 
and -Btag candidates along the boost direction, and the 
known boost of the e + e~ system. A description of the in- 
clusive reconstruction of the B tag vertex is given in [l|| ■ 
For the B° — > K®ir° decay, where no charged particles 
are present at the decay vertex, we identify the vertex of 
the Bqp using the single K® trajectory from the 7r + 7r~ 
momenta and the knowledge of the average interaction 
point (IP) [1(3], which is determined several times per 
hour from the spatial distribution of vertices from two- 
track events. We compute At T and its uncertainty from 
a geometric fit to the T(4S) — * B°B° system that takes 
this IP constraint into account. We further improve the 
sensitivity to At t by constraining the sum of the two B 
decay times (tcp + ttag) to be equal to 2 r with an uncer- 
tainty y2 r, which effectively constrains the two vertices 
to be near the Y(4S) line of flight. We have verified in a 
full detector simulation that this procedure provides an 
unbiased estimate of At. 

The per-event estimate of the uncertainty on At r re- 
flects the strong dependence of the At resolution on the 
K®, flight direction and on the number of SVT layers tra- 
versed by the K° decay daughters. In about 60 % of the 
selected events, each pion track is reconstructed from at 
least one <f> hit and one z hit in the first three layers, 
leading to a sufficient resolution for the time-dependent 
CPV measurement. The average At resolution in these 
events is about 1.0 ps. For events which fail this crite- 
rion or for which At r > 20 ps or the error on At r satisfies 
CAt r > 2.5 ps, the At r information is not used. However, 
since C / can also be extracted from flavor tagging infor- 
mation alone, these events still contribute to the mea- 
surement of Cf and to the signal yield. 

We obtain the probability density function (PDF) for 
the time-dependence of signal decays from the convolu- 
tion of Eq. [1] with a resolution function lZ(6t = At r — 
At, cr At t ) ^ where At is the true value of the proper time 
difference from Monte Carlo. The resolution function is 
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parameterized as the sum of a core and a tail Gaussian, 
each with a width proportional to the reconstructed <7At r , 
and a third Gaussian centered at zero with a fixed width 
of 8 ps [HI . We have verified in simulation that the pa- 
rameters of 7Z(St, crAtJ for B° — > K°tt° decays are simi- 
lar to those obtained from the -Bfl av sample, even though 
the distributions of cta*, differ considerably. Therefore, 
we extract these parameters from a fit to the -Bfl av sam- 
ple. We find that the At T distribution of background 
candidates is well described by a 6 function convolved 
with a resolution function with the same functional form 
as used for signal events. The parameters of the back- 
ground function are determined together with the CPV 
parameters and the signal yield. 

We extract the CPV parameters from an extended un- 
binned maximum-likelihood (ML) fit to kinematic, event 
shape, flavor tag, and decay time variables. We have ver- 
ified that all correlations are negligible, so we construct 
the likelihood from the product of one-dimensional PDFs. 
Residual correlations are taken into account in the sys- 
tematic uncertainty, as explained below. 

The PDFs for signal events are parameterized from a 
large sample of fully-reconstructed B decays in data and 
from simulated events. For background PDFs, we select 
the functional form from from the background-dominated 
sideband regions in our data. 

The likelihood function is defined as: 



C(S f ,C f ,N s ,N B ,fs,fB,a) 



-(N S + N B ) 

AM 



(2) 



II i N sfs4'Ps(x t ,yi;S f ,C f ) + N B f B e c B P B (x l ,y l ;d) 



x [N S (1 - fs)e c sV's{xi;Cf) + N B {l - f B )e%V' B , 

ieb 

where the N selected events are partitioned into two 
subsets: i G g events have At r information, while 
i £ d events do not. fs (/s) is the fraction of signal 
(background) events € g, and the complement to one 
is the fraction of events <= b. The probabilities Vs and 
Vb are products of PDFs for signal (S) and background 
(J5) hypotheses evaluated for the measurements x*i = 

{m B , 7^ miss i 

L 2 /L , cos 6»g, flavor tag, tagging category} 
and iji — {At r , <7At r }- ^'s an< ^ ^'b are the corresponding 
probabilities for events without At r information. In the 
formula, a represents the set of parameters that define 
the shape of the PDFs. Along with the CP asymmetries 
Sf and Cf, the fit extracts the yields Ns and Nb, the 
fraction of events fs and fs, and the parameters a 
which describe the background PDFs. 

Fitting the data sample of 18111 B° -> K°n° 
candidates, we find Ns — 459 ± 29 signal decays with 
S K o w o = 0.40 ± 0.23 and C K o„ = 0.24 ± 0.15, where the 
uncertainties are statistical only. The linear correlation 
coefficient between the CPV parameters is —7.3%. Tak- 
ing into account the selection efficiency and the total 
number of BB pairs in the data sample, we obtain the 
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FIG. 1: Distribution of (a) m m i ss , (b) m B , (c) L^/Lq, (d) 
cos6 B , for background subtracted events on data (dots). The 
solid curve represents the shape of signal PDF, as obtained 
from the ML fit. The insets show the distribution of the data, 
and the PDF, for signal subtracted events. The binning of the 
Z/2 /Lo PDF is coarser where the signal is well separated from 
the background to reduce the number of free parameters. 



branching fraction B(B° -> K° q n°) = (10.4 ± 0.7) x 10~ 6 
which does not include systematic corrections on the 
yield. 

Fig. [T] shows distributions for signal (background) 
events, where background (signal) is subtracted using an 
event weighting technique [19( . Figure shows distribu- 
tions of At r for B and _B°-tagged events, and the asym- 
metry A K o w o(At T ) = [N B o - Ng ] I [N B a + N So ] as a 
function of At r , for background subtracted events. N B o 
(Ng ) represents the number of events tagged as B° {B°). 

In order to validate the IP-constrained vertexing tech- 
nique for CPV measurements we examine B° — * J/ipK® 
decays in data, where J/ifj — > or J/tp — > e + e _ . In 

these events we determine At r in two ways: by fully re- 
constructing the B° decay vertex using the trajectories of 
charged daughters of the J/tp and the Kg mesons (stan- 
dard method), or by neglecting the J/i/j contribution to 
the decay vertex and using the IP constraint and the K® 
trajectory only. This study shows that within statisti- 
cal uncertainties, the IP-constrained Ai r measurement 
is unbiased with respect to the standard technique and 
that the fit values of Sj^ K ° and Cj^ K o are consistent 
between the two methods. 

To compute the systematic error associated with the 
signal yield and CPV parameters, each of the input pa- 
rameters to the likelihood fit is shifted by ±ler from its 
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FIG. 2: Background subtracted distributions of At r for events 
where _B t a g is tagged as (a) B° or (b) B°, and (c) the asym- 
metry A(At r ). The points are data and the curves are the 
PDF projections. 



nominal value and the fit is repeated. Here, ±lcr is the as- 
sociated error, as obtained from the -Bfl av sample (for At 
and tagging) or from Monte Carlo. This contribution to 
the systematic error takes into account the limited statis- 
tics we used to parametrize the shape of the likelihood in 
Eq. [21 We find a systematic error of 0.72 events on the 
yield, and of 0.006 (0.010) on S^o^o (C^o^o). As an ad- 
ditional systematic error associated with the shape of the 
PDF, we also quote the largest deviation observed when 
the parameters of the individual signal PDFs are floated 
in the fit. This gives a systematic error of 11 events on the 
yield, and of 0.019 (0.018) on S K o„ (C K o„ ). The out- 
put values of the PDF parameters are also used to assign 
a systematic error to the selection efficiency of the cuts on 
the likelihood variables. Comparing the efficiency to the 
Monte Carlo simulation we obtain a relative systematic 
error of 3.7%. We evaluate the systematic error coming 
from the neglected correlations among fit variables using 
a set of simulated Monte Carlo experiments, in which we 
embed signal events from a full detector simulation with 
events generated from the background PDFs. Since the 
shifts are small and only marginally significant we use 
the average shift in yield (—2.3 events) and CPV param- 
eters (—0.003 on S K o n o and +0.015 on C K o^o) as the 
associated uncertainty. 

We estimate the background from other B decays to 
be small in the nominal fit. We account for a system- 
atic shift induced on the signal yield and a small system- 
atic uncertainty induced on the CPV parameters by this 
neglected component by embedding simulated B back- 
ground events in the dataset and evaluating the average 



shift in the fit result: +4.5 events on the signal yield, 
+0.003 on S K o v o and -0.002 on C K a„o. We adjust the 
signal yield accordingly and we use half of the shift as a 
systematic uncertainty. 

To quantify possible additional systematic effects, we 
examine large samples of simulated B° —> K®ir° and 
B° — ► J/ipK® decays. We employ the difference in reso- 
lution function parameters extracted from these samples 
to evaluate uncertainties due to the use of the resolution 
function 1Z extracted from the Bft av sample. We also use 
the data-Monte Carlo difference of the resolution function 
in B° — > J/ipK® decays to quantify possible problems in 
the reconstruction of the K® vertex. We obtain a com- 
bined systematic error from this control sample of 0.027 
on S^o^o and 0.006 on C K o^o. 

We assign a systematic uncertainty of 0.002 on S K o^o 
and 0.001 on C K o^o to account for possible misalign- 
ments of the SVT. This does not include the effects as- 
sociated with changes in the At resolution function since 
this is measured on data using the -Bflav control sample. 
We consider large variations of the IP position and reso- 
lution, which produce a systematic uncertainty of 0.004 
on S K o^o and 0.001 on C K o n o. Additional contributions 
come from the error on the known B° lifetime (0.0022 
on both S K o^o and C K o^ ), the value of Amj (0.0017 on 
both S K o^o and C K o v o), and the effect of interference on 
the tag side [H (0.0014 on S K o„o and 0.014 on C K o w o). 
The systematic uncertainties on S^o^o and C K o^ a are 
summarized in Table [U 

For the branching fraction, systematic errors come 
from the knowledge of selection efficiency, (33.6 ± 1.6)%, 
the counting of BB pairs in the data sample, (383.2 ± 
4.2) x 10 6 BB pairs, and the branching fractions of the 
B decay chain B{K% -> ir+ir~) = 0.6920 ± 0.0005 and 
B(n° -> 77) = 0.9880 + 0.0003) 0- The systematic 
uncertainties on the BF are summarized in Table [TXJ 

TABLE I: Summary of contributions to the systematic error 



on S h 



iand Ci. 



AS(1Q- 2 ) AC(10~ 2 ) 



Stat, precision of PDF parameters 


0.6 


1.0 


Shape of signal PDF 


1.9 


1.8 


BB background 


0.3 


0.2 


Correlation among fit observables 


0.3 


1.5 


Vertexing method and !Z(At T , a &t,) 


2.7 


0.6 


SVT alignment 


0.2 


0.1 


Beam spot position calibration 


0.4 


0.1 


B° lifetime 


0.2 


0.2 


Mixing frequency 


0.2 


0.2 


Tag side interference 


0.1 


1.4 


Total 


3.4 


3.0 



In summary, we have performed a measurement of the 
time-dependent CP asymmetries of B° — > K®ir° and 
the branching fraction of B° — > K°ir°. We measured 
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TABLE II: Summary of dominant contributions to the sys- 
tematic error on the branching fraction. 



AB/B(%) 


7r u efficiency 


3.0 


K$ efficiency 


0.6 


Cut on likelihood variables 


3.7 


Shape of signal PDF 


2.5 


BB background 


0.5 


Correlation among fit observables 


0.5 


Number of BB 


1.1 


Resolution function 


0.2 


Total 


5.6 



the CPV parameters C K o w o = 0.24 ± 0.15 ± 0.03 and 
Sk°tt° = 0-40 ± 0.23 ± 0.03, and the branching fraction 
B(B° K°n°) = (10.3 ± 0.7 ± 0.6) x 10~ 6 . The first 
errors are statistical and the second ones are systematic. 
These values are consistent with the Standard Model pre- 
dictions and the experimental value of sin 2/3. The results 
presented in this work supersede previous ones [l3j]. Us- 
ing the rate sum rule from Q and the currently published 
results for the other three B — > Kir modes we find a pre- 
diction for B(B° -> K° s Ti ) sr = (9.0 ± 0.7) x lO" 6 which 
is consistent with our experimental result. Using this re- 
sult the difference between the experimental result and 
the prediction improves from 1.3± 1.1 to 0.9± 1.0, which 
is consistent with zero. 
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